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ABSTRACT 



A highly linear passive floating element wall shear trans- 
ducer has been developed which uses a sensitive frequency modula- 
tion technique to measure a capacitance change proportional to 
the displacement of the floating element for small deflections. 
The floating element is 4 mm square surrounded by a gap nominally 
0.18 mm and protrudes into the boundary layer 9 ym. Scaled on 
inner parameters, each side of the sensing head measures about 
290, the mean gap width is 13, and the element protrusion is 0.7. 
The transducer was tested in a turbulent boundary layer on a flat 
plate at momentum thickness Reynolds numbers of 3140, 3500, and 
3760 and the results compared to wall shear stress estimates ob- 
tained from the law of the wall. The transducer was found to 
predict the steady mean wall shear stress within the limits of 
experimental uncertainty in all cases; transducer measurements 
of mean wall shear stress for momentum thickness Reynolds numbers 
of 3140 and 3500 were within 2 percent of the actual value. 
Transducer measurement of the unsteady mean wall shear stress was 
severely limited by a mechanical resonance in the transducer at 
63 Hz. 



The operating principles and general characteristics of the 
transducer are discussed, its performance is evaluated, and re- 
commendations are made concerning future work. 
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NOMENCLATURE 



surface area of the floating element; 

C capacitance of the cantilever-stationary plate pair 

friction coefficient 

mean spacing between stationary and cantilever plates 
with the cantilever undeflected 



uncertainty or error in the quantity C 

El the flexural rigidity of the cantilever beam 

^ generally, a force; in Appendix C, a functional ex- 

pression for the quantity being measured 

^ frequency in Hz; in Appendix A, a force distributed 

the length of the beam element (f(x,t)) 

center frequency of astable multivibrator 

f j fundamental mechanical resonant frequency of the beam 
element 

G gain in the f requency-to-voltage converter and instru- 

mentation amplifier, V/Hz 

J mass moment of inertia of the floating element 

k surface roughness 

k static spring constant of cantilever beam to transverse 

® loads, 3EI/L ^ 



length of the cantilever plate and beam element 



mass of the floating element (assumed to be a point 
mass) 



mass of the' beam 



the total head, 
dix C) 




also the total error (Appen- 



electrical resistance 



Reynolds number based on x, the distance from the 
plate's leading edge 



momentum thickness Reynolds number 



sensitivity of the transducer in either V/Pa or V/pm; 
also the spacing between capacitor plates (Appendix A) 



thickness of the cantilever plate ; also time 



longitudinal fluctuating fluid velocity 



free stream velocity 



friction velocity 



■ 



transducer offset voltage 



transducer output voltage 



width of the cantilever plate; also (Appendix A) the 
weight per unit length of the beam 



weight of the floating element 
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referring to the flow, the streamwise distance from the 
plate's leading edge; referring to the beam element, 
the distance along the beam from the base of the beam 



referring to the flow, the distance normal to the wall 
measured from the wall; referring to the beam element, 
the deflection of the beam element from its undeflected 
position 



the deflection of the beam with respect to a stationary 
observer 



boundary layer thickness; also (Appendix A) the varia- 
tional operator used in the variational calculus 



•5* displacement thickness 



transverse deflection of the free end of the 
cantilever; also the percent most probable error (Ap- 
pendix C) 



X* 



the viscous length scale. 



V 



U 



beam deflection normalized on mean capacitor spacing d 



normalized coordinate along beam, x/L 



momentum thickness 



transducer inclination angle from the vertical 



floating element angle of rotation 



w 



wall shear stress; mean wall shear stress refers to a 
spatial average as opposed to steady wall shear stress 
which implies a time average 



V 



kinematic viscosity 
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Superscripts 



+ identifies the quantity as being nondimensionalized 

inner coordinates; velocities are normalized on U 
lengths are normalized on A* ^ 



Subscripts 



refers to free stream properties 



w refers to quantities measured at the wall 



in 

and 
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1. INTRODUCTION 



1.1 BACKGROUND 



For many centuries, men have studied the force of drag which 
resists the motion of a body through a fluid. Todd (1967) re- 
ports that as early as the late fifteenth century, Leonardo da 
Vinci conducted tests on three ship models with different hull 
forms, presumably to identify the most efficient design. 

The dawn of the "age of steam” in the late nineteenth 
century motivated researchers to better understand fluid drag in 
order to both predict and minimize its effect. As sailing ships 
were supplanted by larger and faster self-propelled ships, fric- 
tion drag was recognized as a sizeable contributor to the overall 
fluid drag. William Froude (1872) conducted a series of friction 
drag measurements in his towing tank at Torquay, England, to ana- 
lyze the effect of friction drag on ships. Froude built planes 
from wooden planks and towed them edge-first through the water, 
measuring their drag with a dynamometer. In so doing, Froude ac- 
tually determined the mean wall shear stress acting over the sur- 
face of each plane. 

Since 1904, when Prandtl presented his boundary layer theory 
to the Mathematical Congress in Heidelberg, understanding of the 



wall shear stress has grown considerably in sophistication. 
While Froude intuitively grasped that the distribution of wall 
shear stress over a body in a fluid flow was not uniform, theory 
now predicts that the wall shear stress is a function of both 
space and time. This dynamic, nonhomogeneous nature of the wall 
shear stress is particularly evident in the transition and turbu- 
lent regions of the boundary layer. 

The wall shear stress is also an important boundary layer 
parameter in its own right since it is a characteristic of the 
inner region of the boundary layer (near the wall) where viscous 
forces predominate. Many researchers have found that boundary 
layer processes scale on the wall shear stress, particularly 
those related to turbulence in the transition and fully turbulent 
regions of the boundary layer. One example is the familiar "law 
of the wall" in which the turbulent mean velocity profile is non- 
dimensionalized on mean wall shear stress and kinematic viscosity 
[see Coles (1968) and White (1974)]. Burton (1974) found that 
turbulent wall pressure spectra and corresponding fluctuating 
longitudinal velocity spectra scaled on the wall shear stress. 
Gedney (1979) was unable to collapse his transition wall pressure 
spectra because of inaccurate determination of the wall shear 
stress upon which they were nondimensionalized. In short, ongo- 
ing fundamental boundary layer research indicates that more pre- 
cise measures of the wall shear stress than currently available 
are increasingly necessary. 



1.2 WALL SHEAR STRESS 



Wall shear stress, like any stress, is a derived quantity 
and thus cannot be measured directly. The wall shear stress is a 
dynamic quantity defined at a point by the following limit: 

_ lim ^ 

“ A^O A (1.1) 

where F is the friction force acting over area A. Unless the 
measurement of’ the wall shear force is made at a single point, 
the best we can expect to measure is a mean wall shear stress 
averaged over some measurement area. Measurement techniques used 
to determine wall shear stress are, in reality, merely techniques 
which measure quantities related to the wall shear stress from 
which the wall shear stress is inferred. As the measurement area 
approaches a point, the measured wall shear stress approaches the 
true wall shear stress according to equation 1.1. Small measure- 
ment areas thus strengthen the Inference that the measured shear 
stress is the true shear stress, underscoring the need for good 
spatial resolution if the true shear stress at a point is the 
subject of investigation. 
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1.3 MEASUREMENT TECHNIQUES 

Table I, adapted from Winter (1977), provides a brief list- 
ing of the mean wall shear stress measurement techniques current- 
ly in use. No single technique is universally acceptable because 
each has restrictions on its application to different flow condi- 
tions and possesses limited versatility. Further, no technique 
can make an adequate measure of the true wall shear stress, owing 
to excessive spatial averaging or inadequate bandwidth. The hot 
film probe, using the heat transfer analogy, is the only device 
with sufficient bandwidth to measure some of the unsteady char- 
acter of the mean wall shear stress. However, its spatial reso- 
lution and sensitivity are both poor. 



TABLE I 



Mean Wall Shear Stress Measurement Techniques 



Analogies 



Flow About Obstacles 



Heat Transfer 
Mass Transfer 



Preston Tubes 
Stanton Tubes 
Razor Blades 
Steps & Fences 



Liquid Films 



Static Pressure Holes 



Floating Element 



Momentum Balance 



Velocity Profiles (Wall Similarity) 



In an effort to improve the resolution of dynamic wall shear 
stress measurements, a new type of floating element transducer 
has been developed and tested. This device uses a 
capacitance-measuring technique of exceptional sensitivity [first 
developed by Moeller (1983)] to measure the deflection of a small 
floating element in the plane of the wall. The sensitivity of 
the transducer permits its operation with such small motions of 
the floating element that the device exhibits excellent lineari- 
ty. Thus, the powerful techniques of linear systems analysis may 
be used to uniquely and simply relate statistics of the transduc- 
er output to corresponding statistics of the mean wall shear 



stress . 



2 . TRANSDUCER DESIGN 



2.1 GENERAL DESCRIPTION 

The transducer is a passive floating element instrument 
which generates a voltage proportional to the displacement of the 
element and, indirectly, to the shear force acting on the ele- 
ment. This is unlike most floating element devices, which gener- 
ally employ an active force-balance principle so that the element 
remains nearly stationary. The work of Kistler [Paros, (1970)], 
Allen (1980), Frei and Thomann (1980), Escudier, et. al. (1982), 
and Nguyen, Dickenson, et. al. (1984) are all examples of the 
active technique. In the current design, sensitivity has been 
improved to the point that the floating element deflects less 
than 5 pm. For a streamwise deflection of 5 y m, the floating 

element protrudes only 0.63 y m into the boundary layer, which is 
-2 

a y+ of 5x10 . The protrusion due to floating element deflec- 
tion is thus well within the viscous sublayer and can be neglect- 
ed . 



Figure 2.1 presents a schematic of the transducer. The beam 
element is the moving element of the transducer, composed of the 
floating element and a thin cantilever plate. The floating ele- 
ment, mounted in the plane of the wall atop the cantilever plate, 



is sensitive to motion in only one direction. 



Shear forces on 
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the surface of the element in the sensitive direction deflect the 
cantilever plate, producing a change in capacitance with respect 
to a second plate parallel to and mounted a short distance from 
the cantilever plate. This change in capacitance, proportional 
to floating element deflection for small displacements, changes 
the frequency of a square-wave oscillator. The transducer thus 
produces a frequency-modulated signal proportional to the float- 
ing element displacement which is subsequently demodulated and 
filtered by signal conditioning electronics. 



FLOW 










FLOATING 




CANTILEVER 

PLATE- 






FIGURE 2.1 Transducer Schematic - General 



2.1.1 Transducer Mechanical 



The transducer is assembled within a body of 
phenolic-impregnated linen as shown in Figure 2.2. A shallow 
slot, milled in the body halves before assembly, contains a sta- 
tionary capacitor plate of copper-clad laminate epoxied in place. 




FIGURE 2.2 Transducer Schematic 

The movable capacitor plate, made from 0.005 inch (0.13 mm) thick 
sheet brass shim stock, is sandwiched between the halves of the 
transducer body which are epoxied and screwed together. A float- 
ing element 4 mm square is cast in place atop the cantilever 
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plate using an epoxy potting compound. The square shape of the 
floating element was selected to minimize the degradation in 
transducer sensitivity caused by sealing the gap around the 
element . 

The transducer is then mated to an aluminum adapter machined 
to accomodate the test fixture as indicated in Figure 2.3. An 
aluminum cylinder also mates with the adapter, shielding the 
transducer and its electronics from electromagnetic interference. 
Brass shims between the transducer and the adapter are used to 
position the floating element parallel to the top of the adapter 
with the top of the floating element protruding less than 0.1 mm 
from the surface of the adapter. The tops of the adapter and 
floating element, spray painted with flat white enamel before as- 
sembly, were lightly sanded with fine wet sandpaper to reduce the 
misalignment of the floating element with the wall to 9 ym. By 
sanding perpendicular to the transducer's sensitive direction, 
motion of the floating element and curvature of its upper surface 
were minimized. The gap around the element is not symmetric, 
ranging from 0.1 mm to 0.25 mm with a mean of 0.18 mm. An accel- 
erometer mount is provided on the side of the transducer assembly 
inside the cylinder in order to monitor the acceleration of the 
transducer in its sensitive direction. 

Table II summarizes the key lengths of the floating element 
in both dimensional form and in inner coordinates, nondimension- 
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alized on the viscous length A* characteristic of the boundary 
layer in which the transducer was ultimately tested. 




FIGURE 2.3 



Transducer - Pictorial 



w 














TABLE II 

Floating Element Dimensions 



Dimension 




+ 

4 


max. 


0 . 25 mm 


18 


gap min. 


0. 10 mm 


7 


mean 


0. 18 mm 


13 


Protrusion 


9 ym 


0.69 


Element 

Size 


4 mm 


290 



? is based on the experimental result 
U = 1.3 m/s at U = 30.5 m/s. 

T 00 



2.1.2 Transducer Electronics 



The electronics associated with the transducer fall into two 
categories: electronics physically mounted beneath the phenolic 
transducer assembly, and signal conditioning electronics in- 
stalled in a separate enclosure. Each set of electronics uses 
shielded precision power supplies (Analog Devices model 904) 
which are bypassed to ground through 1.0 PF electrolytic capaci- 
tors. Further, each integrated circuit power supply pin is by- 
passed to ground through a 0.1 yF ceramic capacitor. 

Figure 2.4 highlights the significant features of the elec- 
tronic circuitry. 




FIGURE 2.4 Circuitry - Block Diagram 




Transducer-mounted electronics consist of the frequency 
modulator and a buffer amplifier. The frequency modulator is an 
astable multivibrator which produces a square wave nominally 
6.2 V in amplitude at a frequency determined by the variable ca- 
pacitor formed by the fixed and cantilever plates. When the can- 
tilever plate is undeflected, the capacitance measures about 
30 pF and the oscillator's frequency is about 55.6 kHz. An LM310 
voltage follower is used as an isolation amplifier that drives a 
shielded BNC transmission line to the signal conditioning elec- 
tronics. Power is supplied to the transducer-mounted electronics 
through a shielded cable from a power supply in a separate enclo- 
sure. The multivibrator circuit is taken from Jung (1980). 

The signal conditioning electronics are mounted in a shield- 
ed aluminum box containing an integral power supply. An Analog 
Devices 453L precision f requency-to-voltage converter demodulates 
the transducer's frequency modulated signal, converting it to a 
voltage that is proportional to the received frequency. Both the 
gain and the offset of the f requency-to-voltage conversion are 
adjustable for calibration. An AD524 instrumentation amplifier 
(with selectable gains of 1, 10, 100, or 1000) amplifies the de- 
modulated signal and drives both an active filter and a transmis- 
sion line for measuring the fluctuating signal. The active 
filter, taken from Lancaster (1982), is a first-order low pass 
filter with cutoff frequency of 0.2 Hz. The filter drives a sep- 
arate line for measuring the steady transducer output. 



2.2 TRANSDUCER RESPONSE CHARACTERISTICS 



It is shown in Appendix A. 3 that the response of the trans- 
ducer is described by the following equation in the limit of 
small deflections of the floating element: 

3 A G f _ 

V(t) = V + ^ ^ T (t) (2.1) 

8 k d 
s o 

where 

T (t) “ the mean wall shear stress acting on the floating 

w 

element ; 

V(t) - the output voltage from the signal conditioning 
electronics; 

A^ = the surface area of the floating element exposed 

to the flow; 

d^ = the mean separation between the undeflected capa- 

citor plates; 

f^ = the oscillator center frequency; 

G = the total gain in the frequency-to-voltage con- 

verter and the instrumentation amplifier; 

k = the transverse stiffness of the cantilever beam 

s 

when statically loaded at its free end; and 
* the offset voltage. 



Equation 2.1 is applicable for small displacements of the float- 






ing element and for frequencies much smaller than the first reso- 
nant frequency of the beam element. 



Equation 2.1 may be rewritten as 



V(t) = V + s T (t) 

O T W 



( 2 . 2 ) 



where 



s 

T 



9t 



3 A G f 
e o 

8 k d 
s o 



0.40 V/Pa 



(2.3) 



is the transducer's sensitivity. A related sensitivity is the 
transducer output per micron of floating element deflection 



s 

e 




3 G f 

~ 0.56 V/ ym 

8 d 



(2.4) 



The limits of the validity of equation 2.1 may be investi- 
gated by using the predicted beam element deflection shape deter- 
mined in Appendix A. 3 to evaluate the theoretical transducer out- 
put. From the appendix. 



V(t) - 
o 

G f 

o 



1 



0 




1 - 1 
dg 

1 + n(g,t) 



(2.5) 

(A1.15) 



and 



n(C,t) 



( 2 . 6 ) 

(A3. 4) 



2d 



( 3?2 - , 



V(t) - V, 



where 



G f 



■is the nomnalized transducer output andn(C,t)is 



the deflection shape of the beam element normalized on d ^ . 
Equation 2.6 describes the static deflection shape of the beam 
element for a force applied to the floating element and is valid 
for frequencies much smaller than the fundamental mechanical re- 
sonance. The integral in equation 2.5 may be performed numeri- 
cally using equation 2.6 and the normalized transducer output 

plotted as a function of the normalized end deflection of the 

e 

beam -j — . As shown in Figure 2.5, the transducer is linear to 



within 1.25 percent proportional error up to 

2 

that kg * 22.5 N/m, = 16 mm , and d^ - 



= 0.05. Given 



375 ym, the device 



is effectively linear up to 



w 



= 26 Pa since 





k k d 






_ ® 1 1 - so 


e 


T 




— 


W 


A 'max A 


d 




max e e 


o 



max 



(2.7) 



For the turbulent boundary layers considered in the current re- 



search, 



w 



max 



- 5 Pa so that the normal operating regime is the 



much smaller box shown in Figure 2.5. 



The natural frequencies of the transducer response are 
obtained in Appendix B by evaluating the eigenvalues of the beam 
element's free vibration response. For the first two vibration 
modes of the beam element, the rotational inertia of the floating 
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element was calculated to affect the natural frequencies by less 
than 1.0 percent and so was neglected. Calculation of the trans- 
ducer's natural frequencies is complicated by an ambiguity in the 
beam length introduced by the spatial extent of the floating ele- 
ment. Contrary to the model, the floating element is not a point 
mass but rather extends over 6 percent of the cantilever plate's 
length. The fundamental resonant frequency of the beam element 
is estimated from the model to be 61.7 Hz, with an uncertainty of 
6.3 percent. Measurements of the dynamic characteristics of the 
transducer presented in greater detail in section 4.1 indicate 
that the fundamental resonance is actually 63 Hz. 

Characteristics of the transducer are summarized in Table 

III. 



TABLE III 

Transducer Characteristics 



Mass: Floating Element 

Cantilever Plate 



76 mg 
480 mg 



Fundamental Resonance: 



Sensitivity 



Stiffness: 



s 



X 



== 16 g/m 
22.5 N/m 
63 Hz 
0.40 V/Pa 



s 



e 



0.56 V/ym end 



deflection 



Spatial Resolution: 



4 mm X 4 mm 



Inclination of the transducer from the vertical in a gravity 
field generates a response from the transducer, since then a com- 
ponent of the weight of the beam element is transverse to the 
beam. The transducer exhibits buckling behavior when inclined, 
because the transverse force tends to increase the deflection of 
the beam and itself Increases as the local inclination angle of 
the beam Increases. Transducer Inclination response is analyzed 
in Appendix A. 4 in the limit of small static inclinations of the 
transducer base. The expression obtained in the appendix for the 
deflection of the beam element n(^,t)wlth Inclination angle 

0 (equation A4.7) is then inserted in equation 2.5, which is in 

o 

turn integrated numerically to obtain the transducer output as a 
function of inclination angle plotted in Figure 2.6. 

The theoretical response of the transducer to inclination is 
sensitive to the beam length L and mean capacitor spacing d^ , as 
indicated in Figure 2.6 by comparing the solid and dashed lines. 
The solid line is the theoretical transducer response for the 
beam length and mean capacitor spacing used in all calculations; 
the dashed line is the theoretical response for a beam length 

1 mm (3.2%) smaller and mean capacitor spacing 25 ym (6.7%) 
larger. The variations in the parameters used to generate the 
dashed line are representative of uncertainties in those parame- 
ters due to manufacturing limitations. 



A simple experiment was conducted to verify the theoretical 



response predicted in Figure 2.6. Experimental results are also 
plotted in Figure 2.6 and may be seen to lie within the limits of 
uncertainty in predicted response. Also indicated on the figure 
are the expected operating limits of the transducer in the cur- 



rent research from Figure 2.5, corresponding to t. =5.0 Pa at 
zero inclination angle. The transducer output due solely to ro- 
tation of the transducer base one degree from the vertical is of 
the same order as the transducer output due solely to wall shear 
stress, Indicating the importance of maintaining the attitude of 



the transducer. 
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3. EXPERIMENTAL EQUIPMENT AND INSTRUMENTATION 
3.1 EXPERIMENTAL FACILITY 

The transducer was tested in the MIT low-noise, 
low-turbulence open circuit wind tunnel. This facility is shown 
in Figure 3.1 and is described in detail by Hanson (1969). It 
has a flow speed range of 20 to 50 m/s, controlled by varying the 
speed of a DC motor driving the blower. The array of turbulence 
reducing screens between the honeycomb flow straightener and the 
settling zone maintains the free stream turbulence level less 
than approximately 0.05%. A test section 38 cm (15 inches) 
square and 183 cm (6 feet) long was mounted in an open jet confi- 
guration within the airtight blockhouse immediately behind the 
contraction. The test section, lined with smooth sound absorp- 
tive material, has movable walls for controlling the longitudinal 
pressure gradient. A mechanical x-y traverse accurate to 0.025 
mm (0.001 inch) vertically and to 2.5 mm (0.1 inch) horizontally 
is installed on the test section and used to position a hot wire 
or a pressure probe during the measurement of mean flow proper- 



ties in the test section. 



Test Section in 
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3.2 TEST FIXTURE 



The test fixture is a damped aluminium flat plate described 
in detail by Gedney (1979), who used it while investigating tran- 
sition region wall pressure fluctuations. The plate, sketched in 
Figure 3.2, Is 38 cm (15 Inches) wide and 122 cm (48 Inches) 
long; the upper surface is polished to a 0.1 micron root mean 
square (rms) surface roughness. The plate's leading edge is a 
solid tool and jig plate, 1.27 cm (0.50 inch) thick, machined 
into a 6:1 ellipse. The test plate's trailing edge consists of a 
thin splitter plate which extends rearward to prevent coherent 
vortex shedding from creating plate vibrations. 

The test plate was mounted horizontally in the test section 
13 cm (5 Inches) above the bottom to minimize secondary flow ef- 
fects observed by Shapiro (1977). The walls of the test section 
were then adjusted to provide a minimum pressure gradient. The 
variation in pressure obtained over the length of the plate at a 
free stream velocity of 35 m/ s is shown in Figure 3.3. The maxi- 
mum variation observed (excluding the leading edge) was less than 
2 % of the total head. The transducer was installed at the loca- 



tion indicated. 



Splitter Plate 



JO 




FIGURE 3.2 Test Plate Layout 
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Distance from leading edge of plate x (m) 



FIGURE 3.3 Longitudinal Pressure Variation 
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A flow visualization to check the effect of transverse wall 
turbulence was performed using Shapiro's technique, also used by 
Gedney. The results were identical to Gedney's as shown in Fig- 
ure 3.4. The shear stress transducer was tested in the 7.6 cm (3 
inch) hole at the location indicated. 



The boundary layer was tripped to initiate turbulence by two 
thicknesses of duct tape 1 cm in width and 0.5 mm in total 
height. The tape was placed on the upper surface of the plate 
2 cm from and parallel to the leading edge. According to the 
criterion developed by Gibbings [White (1976)], fully effective 
tripping to turbulent flow is assured if 



U k 

oo 

V 



826 



(3.1) 



where k is the height of the surface roughness. In the present 
work, k - 0.5 mm ensures Immediate transition for free stream 
velocity greater than 25 m/s. 
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( 1979 ) ] 



3.3 FLOATING ELEMENT MISALIGNMENT 



Allen (1977) has indicated that the major sources of error 
in floating element shear stress measuring devices are lip forces 
and normal forces which contaminate the friction force that the 
device is being used to measure. Lip forces are developed by a 
nonuniform pressure distribution around the periphery of the 
floating element caused by either flow through the gap around the 
element or misalignment of the floating element with respect to 
the wall. Normal forces introduce error when a nonuniform pres- 
sure distribution over the surface of the floating element shifts 
the normal force's center of effort, generating a moment which 
deflects the floating element. 

Allen found that errors in the measured wall shear stress 
due to lip and normal forces can be larger than the actual wall 
shear stress by as much as a factor of 4. His results indicate 
that proper alignment of the floating element in the wall is the 
key factor in minimizing error, a finding corroborated by Burke 
(1984) who used a floating element to measure skin friction in 
the MIT water tunnel. 

Proper alignment of the floating element in the wall is 
therefore critical in the current design, and was assured by 
measurement with a KD-45 Fotonic Sensor made by Mechanical Tech- 



41 



nology Incorporated. The Fotonlc Sensor is a non-contact measur- 
ing device that projects several beams of light through fiberop- 
tic strands and then measures the quantity of light reflected 
from an object to determine the distance to the object. The Fo- 
tonic Sensor was used in its less sensitive range ( ±2 pm) to de- 
termine that the floating element protruded 9 ym from the wall 
(y^ = 0.69), which was judged to be satisfactory. 



3.4 INCLINATION RESPONSE 



Experimental verification of the transducer's response to 
inclinations from the vertical were made, as shown in Figure 2.6. 
The transducer was attached to the mid-span of a sturdy 1.83 m 
long aluminum I-beam which was raised or lowered at one end a 
known distance to generate inclination angles smaller than one 
degree. An inclination angle of one degree is obtained when one 
end of the I-beam is 32 mm higher than the other end, so it is 
possible to maintain accuracy to about 0.01 degree. The effect 
of I-beam flexure on inclination angle is small to begin with and 
further reduced by placing the transducer at the I-beam's 
mid-span, so that I-beam flexural effects have been neglected. 



w 
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3.5 INSTRUMENTATION 



3.5.1 Dynamic Characteristics 

Transducer dynamic characteristics were investigated using 
the Instrumentation shown in Figure 3.5. The transducer was 
mounted in a phenolic bracket and secured to the vibration ex- 
citer so that the transducer could be shaken either in the plane 
of the wall or normal to the wall. As shown in Appendix A. 2, the 
transducer's response to acceleration of its base and its re- 
sponse to a shear force acting on the floating element are dynam- 
ically equivalent, although caution must be exercised in inter- 
preting results since the acceleration of the transducer's base 
also introduces an effective force along the length of the canti- 
lever plate which is absent during normal operation of the trans- 
ducer. An accelerometer, calibrated using a General Radio model 
1557-A vibration calibrator, was attached to either the phenolic 
bracket or to the transducer body to measure the vibration of the 
transducer body. A Hewlett-Packard 5423A structural dynamics an- 
alyzer was used to generate the random noise excitation of the 
vibration exciter and also to calculate the transfer function 
between the transducer output and transducer body acceleration. 
The spectrum of the random noise excitation was shaped using an 
MXR 31-band equalizer so that the spectrum of the transducer body 
acceleration was approximately uniform over the frequency band of 




FIGURE 3»5 Instrumentation - Dynamic Characteristics 
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interest . 



3.5.2 Test Conditions 

The free stream velocity into the test section was monitored 
by measuring the pressure drop across the wind tunnel contraction 
section with a Betz micromanometer. A calibration was performed 
by Shapiro (1977) to ensure an accuracy of 0.1 m/s. 

The longitudinal pressure variation along the test plate was 
measured using a small static pressure tube (United Sensor and 
Control Corp. PSA~16) and a Validyne DP15~24 differential pres~ 
sure transducer. A Validyne model CD23-1071 electronic manometer 
powered the pressure transducer and provided a digital readout of 
the differential pressure in millimetres of water. 

Mean velocity profile measurements were performed using a 
hot wire anemometer probe (Thermo-Systems Inc. probe model 
1261-T1.5) with a Disa model 55005 battery-powered constant tem- 
perature anemometer and a linearizer built by Moeller using an 
Analog Devices 433 b programmable function generator. The hot 
wire probe was calibrated against a pitot tube used with the Val- 
idyne DP15-24 differential pressure cell. Mean velocity was 
measured on a digital voltmeter after being filtered by a 0.1 Hz 
active low pass filter built by Lueptow (1984). 
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3.5.3 Transducer Calibration 

Wall shear stress transducer electronics were calibrated 
using a General Radio model 1210-C function generator in conjunc- 
tion with a Hewlett-Packard model 5532A electronic counter and a 
digital voltmeter. With gain of the AD524 instrumentation am- 
plifier set at 10, the gain of the frequency to voltage converter 
was adjusted so that a change in input frequency of 500 Hz (from 
55.6 kHz to either 55.1 kHz or 56.1 kHz) produced a change in 
output voltage of approximately 5.0 V. The exceptional sensitiv- 
ity of the transducer made workbench calibration of the entire 
transducer unfeasible, so the transducer was calibrated against 
the law of the wall in situ. 



3.5.4 Wall Shear Stress 

The wall shear stress was measured using the instrumentation 
of Figure 3.6. Coherence between the acceleration of the trans- 
ducer and transducer output was monitored using the structural 
dynamics analyzer to ensure that transducer acceleration did not 
significantly affect transducer output. The structural dynamics 
analyzer was also used in its time record mode to make 
long-period zero drift measurements. 
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FIGURE 3.6 Instrumentation - Wall Shear Stress 
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4. EXPERIMENTAL RESULTS AND DISCUSSION 



4.1 DYNAMIC CHARACTERISTICS 



Transducer dynamic characteristics were evaluated for two 
different transducer configurations: one with a seal in the gap 
surrounding the floating element and another with the gap 
unsealed. For brevity, these two configurations will be referred 
to as the "sealed" and "unsealed" configurations, respectively. 
Investigation of the dynamic characteristics of the two gap con- 
figurations was conducted in two phases. In the first phase, the 
accelerometer used to measure transducer body acceleration was 
attached to the phenolic bracket mounting the transducer to the 
vibration exciter. This test configuration permitted the rota- 
tion of the transducer in the bracket and was used to verify that 
the transducer can discriminate against signals orthogonal to its 
sensitive axis in the plane of the wall, but introduced errors in 
that the accelerometer was not mounted directly on the transducer 
body. The second phase utilized an accelerometer attached di- 
rectly to the transducer body in the location indicated on Figure 
2.2. In this later phase, attention was devoted exclusively to 
improving understanding of the dynamics of the transducer in its 
sensitive direction, the orthogonality of the device having 
already been demonstrated. 
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Phase . The sealed configuration in the initial phase used 
as the gap sealing material an epoxy, ECCOGEL 1365-15, made by 
Emerson and Cuming. Measuring only 15 on the Shore hardness 
scale commonly used in the polymer industry, the epoxy is very 
soft but is nearly incompressible. Owing to the manner of the 
transducer's construction, it was not possible to determine the 
geometry of the gap seal beneath the surfaces of the adapter and 
floating element after the epoxy cured. 

The transfer function between the transducer output and the 
transverse acceleration of the transducer body is shown in Figure 
4.1 (unsealed configuration) and Figure 4.2 (sealed configura- 
tion) . The magnitude plots are normalized on the peak magnitude 
of the unsealed transducer's transfer function in order to high- 
light the order of magnitude loss in response caused by the epoxy 
seal in the gap. The fundamental resonant frequency of the 
unsealed transducer (Figure 4.1) is approximately 63 Hz, very 
close to the 61.7 Hz predicted during design. Figure 4.2 indi- 
cates that the addition of the epoxy seal leads to an increase in 
the fundamental resonant frequency to about 360 Hz as well as an 
apparent reduction in sensitivity. 

The apparent loss of sensitivity is misleading. Table V 
(Appendix B) lists the resonant frequencies for the first three 
vibration modes of the beam element subjected to clamped— free, 
clamped-pinned , and clamped-clamped boundary conditions. The ob- 
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FIGURE 4.1 Initial Transfer Function - Gap Unsealed 
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FIGURE 4.2 Initial Transfer Function 



Gap Sealed 



served fundamental resonant frequency of the sealed transducer 
(360 Hz) is bracketed by that of a clamped~pinned beam element 
(341 Hz) and a clamped-clamped beam element (494 Hz) , neglecting 
the rotational inertia of the floating element which was found to 
be insignificant. This suggests that the beam element in the 
sealed configuration exhibits a pinned condition at the floating 
element or, more correctly, a mixed condition that is a hybrid of 
pinned and clamped conditions but closer to pinned. If such is 
the case, then the sealed transducer does not respond to wall 
shear stress, since the floating element is located at a node of 
the beam element's response and therefore never translates. The 
sealed transducer is thus incapable of responding to forces ap- 
plied to the floating element, but rather responds only to forces 
applied to the cantilever plate. This is shown in Figure 4.2, 

where the transducer is responding to the D'Alembert or inertia 
d^z 



force 




(using the notation of Appendix A. 2) produced on 



the cantilever plate by the transverse acceleration of the trans- 



ducer base. The use of a seal in the gap around the floating 
element is self-defeating in this instance, eliminating transduc- 
er response to the wall shear stress. 



The resonances of both the sealed and unsealed configura- 
tions exhibit considerable damping, approximately 14 percent of 
critical damping in each case. The unsealed configuration's 
response in Figure 4.1a shows an irregularity at 60 Hz caused by 
noise from the field windings of the vibration exciter, while the 



poor response below 10 Hz is due to a low frequency limitation of 
the vibration exciter. The phase angle of both configurations' 
transfer function is anomalous, since it does not change by 180 
degrees as it should when passing through an isolated resonance. 
The irregularity in Figure 4.1a near 50 Hz and the similar irre- 
gularity in Figure 4.2a near 325 Hz are also strange, but may be 
due to the vibration of the bracket securing the transducer to 
the exciter. 



Phase 2 . The second phase of the investigation of transduc- 
er dynamics was undertaken to resolve some of the Issues raised 
concerning the phase of the transfer function and the role of the 
mounting bracket in introducing the anomalies encountered in the 
earlier phase. Also of interest was confirmation of the hypo- 
thesis that the floating element in the sealed configuration was 
effectively held motionless by the seal. In this phase, the ac- 
celerometer used to measure transducer transverse acceleration 
was mounted directly on the phenolic transducer assembly rather 
than on the bracket mounting the transducer to the vibration ex- 
citer. For the ’’sealed” configuration, a single piece of stan- 
dard vinyl electrical tape was securely attached to the surfaces 
of the floating element and adapter to approximate a pinned boun- 
dary condition at the floating element. 



Figures 4.3 and 4.4 illustrate the transfer function 



obtained between the transducer output and the transverse accel- 
eration of the transducer body. The units of the ordinate in the 
magnitude plot are volts per g, where 1 g = 9.8 m/s is the ac- 
celeration due to gravity. 

The unsealed configuration's transfer function (Figure 4.3) 
is virtually identical to Figure 4.1, with the following impor- 
tant exceptions: 

- the phase is well-behaved and very properly shifts 
through 180 degrees when passing through the resonance; 

- the anomaly in the magnitude near 50 Hz has vanished, 
as expected; and 

- the irregularity in magnitude at 60 Hz is sharper than 
that encountered earlier, and gives rise to slight ir- 
regularities at 120 Hz and 180 Hz. 

As discussed earlier, the noise at 60 Hz (and multiples thereof) 
is due to the field windings in the vibration exciter. The 
slight irregularity in magnitude at 100 Hz is due to noise from 
the vibration exciter's blower motor. 

The "sealed” configuration's transfer function is remarkably 
similar to that of Figure 4.2, despite the difference between the 
two seals (electrical tape as opposed to epoxy). In Figure 4.4: 

- the fundamental resonant frequency has increased to 
375 Hz (from 360 Hz in Figure 4.2); 

- the damping has decreased to 7 percent of critical (in 
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FIGURE 4.3 Final Transfer Function - Gap Unsealed 
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FIGURE 4.4 Final Transfer Function - Gap Sealed 



contrast to 14 percent in Figure 4.2); 

- the large irregularity in magnitude just below reso- 
nance has vanished; and 

- the phase is well-behaved and shifts through 180 de- 
grees through the resonance. 

The results of this phase of the investigation indicate that 
the anomalies encountered were indeed experimental and confirm 
the hypothesis that the gap seal immobilizes the floating element 



in translation. 



4.2 WALL SHEAR STRESS 



4.2.1 Law of the Wall 

The law of the wall was used to obtain estimates of the mean 
wall shear stress for comparison with the steady mean wall shear 
stress measured by the transducer. Mean velocity profiles were 
measured above the floating element with a hot wire probe at free 
stream velocities of 26.0 m/s, 30.5 m/s, and 34.5 m/ s correspond- 
ing to momentum thickness Reynolds numbers of 3140, 3500, and 
3760. The velocity profiles were analyzed using the FORTRAN pro- 
gram BLPR0F2 in which the user iteratively guesses the friction 
velocity and the location of the wall (in the measured coor- 
dinates) until the measured data coincides in the inner region of 
the profile with Musker's explicit formula for the turbulent mean 
velocity profile [see Musker (1979)]. BLPR0F2 permits the user 
to fair the measured data into the wall using Musker' s formula 
and then integrates the resulting faired profile using trapezoi- 
dal rule to obtain displacement thickness 6* and momentum thick- 
ness 6 . 



The three mean velocity profiles are plotted in Figure 4.5 
nondimensionalized on inner parameters of the boundary layer. It 
can be seen that the velocity profile for Re = 3500 contains 

U 

little or no wake, whereas the other two profiles exhibit a no- 
ticeable wake. The velocity profile at Re = 3500 was measured 

0 
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FIGURE 4.5 Mean velocity profiles in semilogarithmic coordinates. 

Ordinate scale refers to lowest curve; other curves 
successively displaced upward by 5.0 units. 
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during a different test period than the velocity profiles at 

Re = 3140 and 3760, highlighting the difficulty of reproducing 

9 

test conditions exactly. Since the anomaly in the wake at 
Reg = 3500 does not affect the logarithmic overlap region, the 
change in the outer flow causing the different wake is expected 
to have little effect on the wall shear stress. Because all 
three mean velocity profiles collapse to the classical turbulent 
mean velocity profile, it is assumed that each boundary layer ex~ 
hibits fully developed turbulence despite being tripped into tur- 
bulence so that the application of the law of the wall is valid. 



4.2.2 Transducer Measurement 

Steady . Steady mean wall shear stress was measured by the 
transducer at the same Reynolds numbers as the velocity profiles; 
the results are listed in Table IV, where they are compared to 
the results obtained from the law of the wall. Both the law of 
the wall data and the transducer-measured data are plotted in 
Figure 4.6, taken from skin friction data analyzed by Murlis, 
et. al. , (1982). Transducer measurements at Reg * 3140 and 3500 
are excellent, while that for Re g “ 3760 is a bit high but still 
within experimental uncertainty, estimated in Appendix C to be 66 
percent. The acceleration of the transducer was neglected in an- 
alyzing transducer output, since low coherence between the two 
(at most 10% to 20%) indicated that there was little systematic 



relationship between the transducer output and acceleration of 
the transducer body. 



TABLE IV 
Results 





3140 


350C 


3760 


T 1 

law of the wall 


1.25 Pa 


2.05 Pa 


2.20 Pa 


^measured 


1.3 Pa 


1.7 Pa 


3.1 Pa 



Transducer zero drift was surprisingly low despite the 
transducer's extreme sensitivity, averaging 20>4 of the mean sig 
nal over the period of the measurement (1 to 1.5 minutes). Zero 
drift was measured by covering the floating element with an index 
card while the wind tunnel was running at the proper speed. A 
short length of wire 0.5 mm in diameter was taped beneath the 
card as a spacer to keep the card from contacting the floating 
element. Sample plots of the zero drift measured by the struc- 
tural dynamics analyzer over a period of 1 minute are shown in 
Figure 4.7 where they are compared to the mean output voltage de- 
termined by subtracting the average reference voltage from the 
transducer total output voltage. Each data run took 3 to 4 
minutes from the start of the initial zero measurement to the 
measurement of the wall shear stress and the concludxng zero 



measurement . 
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FIGURE 4.6. Skin friction coefficient variation with momentum thick- 
ness Reynolds number [^from Murlis, et. al. , (1982)J : +, 
measured data using wall shear transducer; o, data from 
law of the wall. 
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FIGURE 4. 7 Zero Drift 
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Unsteady . The unsteady mean wall shear stress spectrum for 
Reg = 3140 and 3760 is shown with the background spectrum in Fig“ 
ure 4.8. Although the plots show some dynamic character of the 
unsteady wall shear stress, they show little of its structure 
owing to the low resonant frequency of the transducer. If it is 
assumed that an instrument is usable up to a frequency one order 
of magnitude less than its first natural frequency, then the 
transducer can only be used up to 6.3 Hz. Burton (1974) indi- 
cates that the frequency related to sublayer bursting in a turbu- 
lent boundary layer is given approximately by 

= 0.03 (^-1) 

oo 

where 6* is the displacement thickness. For the conditions under 
which the transducer was tested, equation 4.1 predicts the burst 
frequency to be about 300 to 500 Hz, well above the nominal 



bandwidth of the transducer and thus unmeasurable. 
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FIGURE 4.8 



Wall Shear Stress Spectral Density 
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5. RECOMMENDATIONS 



In the short term, the wall shear transducer is usable for 
measuring the steady mean wall shear stress. The sensitivity of 
the device is so great, however, that calibration must be per- 
formed ^ situ , preferably in a laminar boundary layer where the 
wall shear stress is well known. The electronics should be up- 
graded by the installation of precision instrumentation 
amplifiers to replace the commercial-grade amplifiers currently 
in use. This should reduce the high frequency noise observed at 
the transducer output and improve the stability of the astable 
multivibrator, whose center frequency was observed to drift from 
55.6 kHz to 46.5 kHz over seven days of continuous operation. 



Long range plans for the transducer should include increas- 
ing its bandwidth, decreasing its sensitivity to inclination, and 
gaining an understanding of how the gap around the floating ele- 
ment affects the shear stress measurement (given that the gap 
cannot be sealed, as the current research suggests). 

As regards the bandwidth, the first resonant frequency of 
the beam element needs to be Increased at least two orders of 
magnitude so that the transducer can measure frequencies near the 
burst rate. A larger bandwidth is desirable so that the device 
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can always be used more than an order of magnitude below its fun- 
damental resonance, which is far less complicated than attempting 
to use an instrument through its first resonance. The transducer 
should be as simple to use as possible. The linearity of the 
transducer, coupled with sufficient bandwidth, will make it a 
powerful instrument for measuring the dynamics of the wall shear 
stress. As linear system theory demonstrates, only in a linear 
system can the statistics of an output (such as transducer output 
voltage) be related definitively and unambiguously with statis- 
tics of the input (such as wall shear stress). 

Next, the sensitivity of the transducer to inclinations from 
the vertical should be reduced by at least two orders of magni- 
tude so that inclination has a second-order effect on transducer 
output (it is assumed that the steady transducer signal is 
zeroth-order while the unsteady signal is first-order) . Coupled 
with the need to maintain roughly the same sensitivity to wall 
shear stress s^, the requirements for an improved physical design 
imply that the mass of the floating element and capacitor plate 
must be reduced several orders of magnitude. If the size of the 

floating element is reduced, then equation 2.3 indicates that 
f G 

Q must be increased a compensating amount to maintain the 
^s ^o 

desired s^. Decreases in require further decreases in mass in 
order to maintain the required bandwidth. Increases in G lead to 
only marginal improvements in transducer performance since noise 
is generally amplified along with the signal. 
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Finally, and most important, a better understanding of the 
role of the floating element gap in affecting the boundary layer 
must be reached. There is little in the literature that sheds 
any light on this issue. Allen (1977) identified that small gaps 
Increase the sensitivity of the floating element transducer to 
errors introduced by misalignment of the floating element with 
the wall. Winter (1977) mentions that the effect of gaps perpen~ 
dlcular to the flow is poorly understood, but that the effect of 
the gap should be negligible if the gap width is smaller than 100 
in wall coordinates. McMichael, et. al. , (1979) found that a 
transverse slot in the wall of a turbulent flow affects the boun~ 
dary layer downstream of the slot. In the authors' words, "Evi- 
dently, temporary release of the no slip condition as the bounda- 
ry layer passes over the slot results in a local decrease in tur- 
bulence production as indicated by the weakened intensity near 
the wall just downstream of the slot." (p. 417). Since the gap 
is adjacent to the wall area where the shear stress will be meas- 
ured, it is unclear whether or not the boundary layer will "heal" 
itself sufficiently to permit an unbiased measure of the wall 
shear stress, particularly if the measurement area is small. The 
effect of the gap is a fundamental issue which has not yet been 
adequately addressed. 



The most fruitful course of action appears to be the micro- 
miniaturization of the transducer to a scale where a very small 
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gap and precise alignment of the floating element can be mainta- 
ined. Peterson (1982) discusses the merits of unicrystalline 
silicon dioxide at great length; small cantilever beams 8 ym 
long, 2 ym wide, and 0.1 ym thick with a fundamental resonant 
frequency of 1.2 MHz have already been made from the material. A 
microminiaturized wall shear stress transducer with such a small 
cantilever, using the capcitance-measuring technique outlined in 
the current research, has a number of distinct advantages. A 
floating element with sensing area 2 ym by 0.1 ym would have su- 
perb spatial resolution, excellent sensitivity, and a large 
bandwidth. The considerable reduction in mass would greatly re- 
duce the sensitivity of the device to inclination. Further, it 
should be possible to keep the gap around the floating element on 
the order of about 10 angstroms so that the gap has a negligible 
effect on the flow. Precise alignment of the floating element is 
possible on the microscopic scale, which may circumvent the ex- 
treme sensitivity of the floating element to misalignment-induced 
errors when the gap is small. Recent trends suggest that such a 
transducer could be mounted on a 20 ym chip with its own acceler- 
ometer and compensation network, permitting an array of such 
transducers to measure the spatial distribution of the shear 
stress over the surface of the wall. 
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6. CONCLUSIONS 



The development of a passive floating element wall shear 
stress transducer , using a sensitive capacitance~measuring 
technique has been presented. In particular, the viability of 
using a frequency modulation scheme to perform the capacitance 
measurement has been demonstrated. The electronic sensitivity is 
such that the transducer is able to measure a 0.1 ym displacement 
of the floating element. Expected motions of the floating ele- 
ment are less than 5 ym and should not adversely affect the char- 
acter of the boundary layer in which the measurement is made. 
The device is highly linear, enhancing its attractiveness for dy- 
namic wall shear stress measurements. 



The 


transducer has been 


tested 


and 


found 


to predict 


the 


steady 


mean wall shear 


stress 


to 


within 


the limits 


of 



experimental uncertainty (66%); two of the three data points ob- 
tained were within 2 percent of the actual steady mean wall shear 
stress. The transducer must be calibrated in situ in order to 
reduce the experimental uncertainty; workbench calibration is 
not feasible because of the transducer's exceptional sensitivity. 
The current transducer is not practical for dynamic measurements 
of the mean wall shear stress since it does not possess adequate 



bandwidth. 



The gap around the floating element has been verified as 
being the greatest single vulnerability of the floating element 
transducer. It is not possible to seal the gap, since a gap seal 
in the current design immobilizes the floating element and pre- 
vents transducer response to wall shear stress. On the other 
hand, an open gap affects the boundary layer and Introduces er- 
rors of its own to the wall shear stress measurement. 
Fundamental research is needed in this area in order to fully 
understand the gap effect. 



Finally, a course of action has been recommended to improve 
the capabilities of the wall shear transducer so that it is prac- 
tical for dynamic measures of the shear stress. In order to in- 
crease bandwidth, reduce transducer sensitivity to inclination, 
and minimize the adverse effects of the gap, the transducer 
should be microminiaturized. Microelectronics technology is be- 
coming increasingly capable, and can already maintain the neces- 
sary sub-micron tolerances. 
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APPENDIX A 
CALCULATIONS 



A.l Transducer Model 

A. 2 Beam Element Equation of Motion 

A. 3 Transducer Response 

A. 4 Transducer Sensitivity to Inclination 



A.l Transducer Model 



Using the model of Figure A.l, it is possible to develop a 
general theory to describe the output of the wall shear transduc- 
er with arbitrary input. Figure A.l highlights the important fa- 
cets of transducer operation: 

- the: wall shear stress causes a deflection of the beam 
element; 

- deflection of the beam element alters its capacitance 
with respect to a stationary parallel plate; 




FIGURE A. 1 Transducer Model 




the change in capacitance 



of 



the 



beam 



element-stationary plate changes the oscillation fre- 
quency f of an as table multivibrator, since the oscil- 
lation frequency is inversely proportional to the capa- 
citance; and 

the frequency difference between the oscillator signal 

and a reference frequency f is then converted to a 

ref 

voltage which is added to an offset voltage V , gener- 



ating the transducer output. 



The capacitance of a parallel plate capacitor is given by 



where K is a constant of proportionality, L is the length and w 
the width of each plate, and s is the distance between plates. 
If one plate is deflected a small transverse amount, then the 
distance s is no longer constant but rather a function of x, the 
coordinate along the undeflected length of the movable plate 
measured from the wall. Further, the deflection may also be a 
function of time. Effectively, s not a constant leads to an in- 
finite number of parallel capacitances 



o 



s 




(Al.l) 




(A1.2) 



in the limit of small plate deflections. If now s(x,t) is de- 



fined as 



s(x,t) = + y(x,t) 



(A1.3) 



where d^is the undeflected capacitor plate spacing and y(x,t) is 
the transverse deflection of the movable plate, then 



dC = K 



w 

d + y(x,t) 
o 




(A1.4) 



Nondimensionalizing and defining 



X 

L 



= C 



^ = n(C,t) 



(A1.5) 



and noting that 




(A1.6) 



where is the capacitance of the stationary-movable plate pair 
with the movable plate undeflected, equation Al.4 becomes 



dC 



C 

o 



dC 

1 + n(C,t) 



(A1.7) 



Equation A1.7 is valid only for small deflections n(5,t) , for 
which three dimensional effects are insignificant. The total ca- 
pacitance may be found to be 
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/ 



C = / dC = C 



AL 

o / 1 + n(?,t) ’ 
0 



(A1.8) 



’summing" over all the infinltesmal capacitances dC. 



The oscillation frequency of the astable multivibrator is 
given by 



f 



1 

R C 



(A1.9) 



where C is the capacitance of the transducer and R is a feedback 
resistance (constant) . Then, for a deflection of the movable 
plate n( t), the oscillator frequency is 



f 



1 

R C 



RC 



1 



0 




dC 

1 + n(Cjt) 



(Al.lO) 



or 



f = 



f 

o 



1 




d^ 

1 + n(?,t) 



0 



(Al.ll) 



Since 



V(t) = V + G 
o 



( f - f 



ref 



it is appropriate to select f ^ = f to obtain 

ref o 



(Al. 12) 



I 







V(t) = 

o 



+ G f 

o 




- 1 



(A1.13) 



or 



V(t) = V + G f 
o o 



1 




1 + n(C,t) 



(A1.14) 



using equation Al.ll. 



The normalized transducer output is given by 



V(t) - 



G f 



J 1 + n(C, 



- 1 



(A1.15) 



t) 



0 



The deflection n( C,t)is determined by the transducer's excitation 
and dynamics considered in the following section. 



A. 2 Equation of Motion 



The dynamics of the movable plate of the transducer deter- 
mines how the transducer will respond to applied loads. The mov- 
able plate Is cantilevered at Its base and has a floating element 
mounted on Its free end. The movable plate and floating element, 
together called the beam element, may be modeled as In Figure 
A. 2. 




F(t) 



M, J 



y(x,t) = z(x,t) - z(0,t) 



z(x,t) 



z(0,t) = z (t) 
o 




FIGURE A. 2 Beam Element Dynamic Model 



oo 

<~>s: 
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The floating element is modeled as a mass M and mass moment 
of inertia J lumped at the free end of the cantilever plate, mo- 
deled as a uniform beam of flexural rigidity El and mass per unit 
length y. The motion of the beam element is described by z(x,t) 
in a stationary reference frame. In order to develop the equa- 
tions of motion for general excitation, a concentrated transverse 
force F(t) is applied to the floating element and a distributed 
transverse force f(x,t) is applied along the beam. 



Of interest in the current situation is the motion of the 
beam element relative to its base rather than its absolute motion 
z(x,t), since the stationary capacitor plate is actually fixed 
with respect to the cantilever's base. In effect, it is desired 
to obtain the equations of beam element motion in the 
non-inertial reference frame translating with the cantilever's 
base. 



If 



z^(t) = z(0,t) 



(A2.1) 



is the translation of the base of the beam element , then the 
relative motion of the beam with respect to its base is 



w 



'■ ^ 






# ^ 
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y(x,t) = z(x,t) - z^Ct) . 



(A2.2) 



Using the variational technique [Crandall, et. al. (1968)] and 
assuming small deflections of the beam, the equation of motion of 
the beam element is obtained as 



M 




d^z 



dt" 



f (x,t) 



(A2.3) 



subject to the boundary conditions 




(A2 ,4a) 
(A2.4b) 

(A2.4c) 
(A2. 4d) 



where 6 is the variational operator. For the case where the beam 
is cantilevered at its base. 



6y 



x=0 




(A2.5) 



so that boundary conditions A2.4a and b are identically satis- 
fied . 



•iRlfA 




Clamp ed- Free Beam. 

In the case of normal transducer operation, the beam element 

3y 

is free at x=L so that the variations in y and at x=L do not 
vanish identically, requiring the terms in brackets in equations 
A2.4c and d to be equal to zero. Thus, 



a2y ; .. aS 






dt^ 



(A2 . 6) 



with boundary conditions 



y(0,t) = 0 



iz 

3x 



(0,t) E 0 



(A2.7a) 

(A2.7b) 






M 



d^ z 

dF" 



(A2.7c) 



FT 



(A2.7d) 



Note that acceleration of the transducer base introduces 



D'Alembert or inertial force excitations proportional to 



o in the field equation (A2.6) and the force boundary condi- 

I 

tion (A2.7 c). If the transducer base is non-accelerating 
d^z 

( = 0 )> the transducer is kept vertical so that f(x,t) = 0, 

dF 

and the floating element is properly aligned so that there are no 

error forces on the floating element (F(t) = t (t)A ) , then the 

w e 



field equation becomes 



ou 



3^v 



(A2.8) 



subject to the revised boundary condition (A2.7c) 



= 7(t)A 

8t^ w 



(A2.9) 



with the other three boundary conditions unchanged. Therefore, 

in the absence of the error terms specified above, the wall shear 

stress excites the homogeneous response of the beam y(x,t) 

through the force boundary condition, equation A2.9. Transverse 

acceleration of the transducer base introduces not only a 

d^ z 

D'Alembert force at the floating element (x=L) equal to -M — 7 -° > 

d^z 

but also a distributed D'Alembert force -y along the length 

dt 

of the beam. Thus, the response of the beam element to 
transverse acceleration of its base is dynamically analogous to 
its response to the wall shear stress except that transverse ac- 
celeration gives rise to the aforementioned distributed 
D'Alembert force. 



Clamped - Plnned Beam . 



If the beam element is pinned at the floating element, then 
= 0 at x=L and the equations describing the motion of the beam 
element are 



P 




f (x,t) 



d^z 

o 



(A2.6) 



with boundary conditions 



y(o,t) = 0 

9x 



(A2.7a) 

(A2.7b) 



y(L,t) = 0 



(A2. 10) 






. aV(L,t) 



(A2.7d) 



Because F(t) does not excite the transducer, the wall shear 
stress cannot either. Thus, if the floating element exhibits a 
pinned end condition, the beam element will not respond to the 
wall shear stress although it will respond to the tranverse ac- 
celeration of the transducer base (equation A2.6 applies). 
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Clam ped- Clam ped Beam . 

If the floating element exhibits a clamped end condition, 
then equation A2.6 is still the appropriate field equation, with 
y and identically zero at x=0 and x=L. Again, the beam ele- 
ment will not respond to the wall shear stress but will respond 
to transverse acceleration of the transducer base. 
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A. 3 Transducer Response 



For frequencies much smaller than the fundamental resonant 
frequency of the beam element, the deflection shape of the beam 
element is satisfactorily described by its static deflection 
shape. Assuming that the transducer base is non-accelerating, 
that the transducer is kept vertical, and that there are no error 
forces on the floating element, equation A2.8 may be solved with 
equations A2.7a,b, and c and equation A2.9 for the static deflec- 
tion shape (which is not a function of time). Assuming 

F(t) = A T (t) (A3.1) 
e w 



leads to 

A _ 

(3^2 _ ^3) (A3.2) 

s o 



where k = 3 EI/l31s the effective stiffness of the beam in re- 
s 

sisting an applied transverse force at the free end. 
alternatively, since 

A _ 

T^(t) • (A3.3) 

s 



is the transverse deflection of the free end of the beam. 
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n( 5 >t) - ( 3{2 - 53 ) j(|.) 



describes the deflection shape of the beam element. 



V(t) - V 



G f 



f- 

/ 1 + n(C,t) 



- 1 . 



0 



it is possible to evaluate the transducer output as 






A. 3.1 Linearized Response 



If 



n(C,t) 



« 1 , then 



1 



/ 



1 , 1 

- = 1 + 

1 - I n(£.t) dC 0 



dC 



1 + n(C,t) 



■ flu,< 
0 



Thus , 



V(t) - V 



G f 



fnU, 



t) dC 



from which, with equation A3.2, 



V(t) - V A 

o e 



G f 



T (t) r 

/ 

: d 



- C^)dC 



2 k d 

s o 0 



(A3. 4) 

Since 

(Al. 15) 



a function of 



dC . (A3. 5) 



(A3. 6) 



(A3. 7) 



Upon integration, equation A3. 7 leads to 



(A3. 8) 



yi 



v(t) - V, 



G f 



3 A 



2 k d 
s o 



— T (t) 
W 



which is shown in Figure 2.5 as a function of 



A T 
e w 

k d 
s o 



(A3. 9) 



the normalized end deflection of the beam. 



A. 3. 2 Exact Response 

Using equation A3. 4, the integral in equation Al.15 can be 

e 

computed numerically for different . The results of this in~ 

o 

tegration are shown in Figure 2.5. 
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A. 4 Sensitivity to Inclination 



When the transducer is inclined from the vertical by an 
angle 9^ , the beam element deflects by its own weight. For 
small inclinations from the vertical, the transverse component 




FIGURE A. 3 Transducer Inclination Geometry 




f^’‘v 




of a load F is given by 



F^ = Fsini|;-F4;, 



(A4.1) 



where 'P is the local angle between the tangent to the beam and 
the vertical shown in Figure A. 3. Now, 



where 6 = is the local angle between the tangent to the beam 
and the undeflected beam position and 6^ is the angle the tran- 
ducer base is inclined from the vertical. 

The beam element consists of a floating element of mass M 
and a cantilever plate of mass per unit length y . In a gravity 
field, these possess weights of W and w, respectively. If the 
transducer is inclined from the vertical with no shear stress ap- 
plied to the floating element and without accelerating the base 
of the transducer transversely, then, from equations A2.6 and 
A2.7 and by considering only static deflections y(x) of the beam 
element , 



ip = d + e 



0 




(A4.2) 



El - f(x,t) 



(A4.3) 



subject to the boundary conditions 



(A4.4a) 



y4 



y(0) E 0 

iz(o) = 0 

dx (A4.4b) 

-El ^ = F(t) (A 4 . 4 c) 

d^y (L) = Q ^ (A4.4d) 

dx^ 

Both f(x,t) and F(t) are developed by the weight of the canti- 
lever plate and the floating element according to equations A4.1 
and A4 . 2 : 



f(x,t) - w ( 9 + ) (A4.5) 

0 dx 

F(t) = w ( + ^ ). (A4.6) 



Equations A4.3 through A4.6 are then solved to obtain 



n(C) ~]Ci + C2 exp (Be) + 



exp(- -| O 





. . 1 


1 ) e L 

f n 


C3 sin 2 


BC + cos -y BC 


> 0 

1 



(A4.7) 



wLi ^ 

where 6 is in radians and B^ = ~ — • The coefficients C were 
0 El i 

calculated using the FORTRAN program GRAV3 , owing to the complex- 
ity of the algebra and in order to facilitate the ready calcula- 
tion of deflection shapes for different transducer parameters. 
The transducer parameters used for the calculations are: 

L = 31 nnn, 
do = 375 ym , 
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M = 76 mg, 

= 16 g/m, and 
El = 2.84x10 Nm^ 



for which 

= 1.1568 

C = 0.9807 

2 

C = 2.1674 

3 

C = -2.1375 

4 

and 3 = 0.254314. 



APPENDIX B 



TRANSDUCER NATURAL FREQUENCIES 



B . 1 Clamped - Free Beam 

The homogeneous form of equations A2.6 and A2.7 may be 
solved in order to determine the natural vibration frequencies of 
the clamped-free beam element: 



y 




(Bl.l) 



subject to the boundary conditions 

y(0,t) = 0 (B1.2a) 

s 0 (B1.2b) 

3x 

Separating variables and assuming time harmonic behavior leads to 

y(x,t) = X(x) T(t) (B1.3) 



where 
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X(x) - Cl sinh Xx + C 2 cosh Xx + C 3 sin Xx + cos Xx (Bl.4a) 

and T(t) = (Bl.4b) 

with 

X^^ = ^. (B1.5) 

El 



Imposing the boundary conditions (Bl.2) obtains 



cosh XL + cos XL - y(sinh XL - sin XL) ^ M 
sinh XL - sin XL - y(cosh XL - cos XL) y 



(Bl. 6 a) 



where 



Y 



sinh XL + sin XL 



cosh XL + cos XL 



0 )^ J 
XEI 



g)^J 

XEI 



(cosh XL 



(sinh XL 



- cos XL) 



+ sin XL) 



(Bl.6b) 



The roots of equations Bl .6 determine the eigenvalues X^ 

The first three natural frequencies have been calculated and are 
tabulated in Table V. The effect of the floating element's mass 
moment of Inertia was approximated by modeling the element as a 
slab 4 mm by 4 mm and 2 mm thick, and was found to be negligible 
for the first three vibration modes. 
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B . 2 . Clamped - Plnned Beam 

The determination of the resonant frequencies of the 
clamped-pinned beam follows the procedure for the clamped-free 
beam except that equation Bl.2c is replaced by 

y(L,t) E 0 . (B2.1) 

The characteristic equation then becomes 

sinh XL - sin XL - yCcosh XL - cos XL) = 0 • (B2.2) 

where y is as given in equation Bl.6b. The natural frequencies 
of the first three clamped-pinned modes is tabulated in Table V, 
both including and neglecting the effect of the rotational iner- 
tia of the floating element. Blevins' (1979) tabulation of the 
natural frequencies for the clamped- pinned beam are the same as 
the natural frequencies predicted from equation B2.2 above. 



B.3 Clamped" * Clamped Beam 



For the clamped-clamped case, the mass moment of inertia no 
longer enters the problem since 

_^(L,t) i 0 (B3.1) 

3x 

is the new boundary condition that replaces equation Bl.2d. The 
natural frequencies of the first three clamped-clamped modes are 
taken from Blevins (1979) and shown in Table V. 



TABLE V 



Summary of Beam Element Natural Frequencies 



Clamped- Free Beam 


b 


b 


b 


J = 0 


61.7 


414.2 


1202 


J ^ 0 


61.7 


411.2 


1180 


A% 


0 


0.7% 


1.9% 


Clamped-Pinned Beam 


b 


b 


b 


J = 0 


340.6 


1104 


2303 


J 0 


339.4 


1091 


2246 


A% 


0.4% 


1.1% 


2.5% 


Clamped-Clamped Beam 


b 


b 


b 




494.2 


1362 


2671 



All frequencies are in Hz. 

-4 2 

J 7^ 0 indicates J assumed to 1.267x10 mg m , 
modeling the floating element as a rectangular slab 
4 mm X 4 mm X 2 mm. 

The floating element mass M and (where appropriate) 
mass moment of inertia J are modeled as being lumped 
at the undamped end of the beam element. 

Clamped-clamped data and verification of clamped-pinned 
calculations are from Blevins (1979). 



Notes: 1. 

2 . 

3. 

4. 
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APPENDIX C 



UNCERTAINTY ANALYSIS 



From equation 2.1, 



8 k d 

X = s_o_ j 

3 A f G ° 

e o 



(C.l) 



It Is shown in Cook and Rabinowicz (1963) that if 



q = F(a.) 



(C.2) 



where F is a function of n independent variables a^ , then the 

2 

most probable mean square error e^ in the quantity q is 



■ t (-f )■ -I . 

i=i ^ ^ 



(C. 3) 



assuming the a^ are statistically independent. Each term 
e^ denotes the error or uncertainty in the quantity C . From 
equation C.3, 



V 






0\ T 

W T, 



(C.4) 



e . 

where e.E — and the uncertainty in V(t) has been neglected since 
"" ^i 

it is an order of magnitude smaller than e 

Vo 



3wt ~ 

Since k = 3EI/L^ = j and A = w^ , equation C.3 obtains 



4L 



iUZ 



€ 



2 





(3C^)^ 






(C.5a) 



(.05)^ + (.05)^ + 9(.01)^ + 9(.05)^ =.0284 (C.Sb) 



and 



= 2 (.05) 2 = 0.0050 



(C.6) 



so that, from equation C.4, 



0.0284 + (.05)^ + 0.005 + (.01)^ + (.05)^ + 

(.10)^1 h-LX \2 =. 0.429 

\(1 Pa) (.4 V/Pa)/ 



(C.7) 



and 



e ~ 

T 






0.66 



(C.8) 



\nh'2. ^ I 
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